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1 REFACE 


This  report  describes  the  work  performed  under  the  DOT/FAA  Core  Engine 
Noise  Control  Program  (Contract  D0T-FA72WA-3023) . The  original  work  under 
this  contract  is  in  Report  Number  FAA-RD- 74-125,  Volumes  I, II,  and  III. 

* 

Supplements  to  Volumes  II  and  III  report  additional  work  undertaken  under 
this  program  after  completion  of  the  work  reported  in  the  original  three 
volumes. 


The  objectives  of  the  program  were: 

• Identification  of  component  noise  sources  of  core  engine  noise 
(Phase  I). 

• Identification  of  mechanisms  associated  with  core  engine  noise 
generation  and  noise  reduction  (Phases  II  and  III) . 

• Development  of  techniques  for  predicting  core  engine  noise  in 
advanced  systems  for  future  technology  aircraft  (Phase  IV) . 

• Extension  of  the  core  noise  prediction  (Phase  V) . 

• Update  of  ; he  core  engine  noise  control  effort  (Phase  VI) . 


The  objectives  were  accomplished  in  six  phases  as  follows: 

• Phase  I - Analysis  of  engine  and  component  acoustic  data 

to  identify  potential  sources  of  core  engine 
noise  and  classification  of  the  sources  Into 
major  and  minor  categories. 

• Phase  II  - Identification  of  the  noise  generating  mechanism 

associated  with  each  source  through  a balanced 
program  of: 

Analytical  studies 

Component  and  model  tests 

Acoustic  evaluation  of  data  from  existing 
and  advanced  engine  systems. 


• Phase  III  - Identification  of  noise  reduction  mechanisms  for 

each  source  through  a program  with  elements 
similar  to  Phase  II. 


• Phase  IV  - Development  of  improved  prediction  techniques 

incorporating  the  results  obtained  during  the 
preceding  two  phases. 


ill 


• Phase  V - Analysis  of  low  frequency  noise  transmission 

through  turbine  blade  rows  and  addition  of  engine 
and  component  data  to  the  prediction  method  for 
core  noise. 

• Phase  VI  - Analytical  studies  of  turbine  source  noise  sup- 

pression and  parametric  trends  in  turbine  tone/ 
jet  stream  interaction;  an  experimental  study  of 
compact  low  frequency  noise  suppressors  and  a 
prediction  model  update. 


The  work  accomplished  is  reported  in  five  volumes  correspond ing  respec- 
tively to  the  five  objectives  stated  above. 


Volume  I 
Volume  II 

Volume  III 

Volume  III 
Supplement  I 


- Identification  of  Component  Noise  Sources 
(FAA-RD-74-125,  I). 

- Identification  of  Noise  Generation  and  Suppression 
Mechanisms  (FAA-RD-74-125,  II). 

- Prediction  Methods  (FAA-RD-125,  III). 

- Extension  of  Prediction  Methods.  (FAA-RD-74-125,  III-I) 


• Volume  II 

Supplement  I - Extension  to  Identification  of  Noise  Generation  and 
Suppression  Mechanisms.  (FAA-RD-74-125,  II-I) 


A visual  representation  of  the  overall  program  and  report  organization 
is  shown  on  page  v and  vl. 

This  volume  documents  the  results  of  the  Phase  V activity  which  extended 
the  core  noise  prediction  technique  through: 

1)  Examination  of  additional  engine  data  to  validate  the  prediction 
method. 

2)  Re-evaluation  of  the  prediction  at  high  power  settings  through 
excmination  of  engine  data. 

3)  Analytical  determination  of  turbine  low  frequency  noise  attenuation 
in  order  to  add  component  data  to  the  engine  prediction  lines. 
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SUMMARY 


The  DOT/PAA  Core  Noise  Prediction  1 d^l  was  modified  to  inciud'  the 
results  of  additional  engine  noise  data.  The  three  line  prediction  method  was 
improved  to  a single  prediction  line  method  by  utilizing  a turbine  work  extrac- 
tion factor.  The  improvement  in  the  prediction  model  has  yielded  a correla- 
tion within  ±2  dB  for  all  data  analyzed  to  date. 

A total  of  thirteen  sets  of  engine  data  were  evaluated.  Three  of  these 
sets  were  from  non-General  Electric  teste:  a JT8D-109  (P&WA)  low  bypass  ratio 

turbofan,  a TPE  331-201  (Garrett)  turboshaft  engine  and  a TF34  (NASA)  high 
bypass  turbofan.  The  spectra  and  directivity  of  all  three  sets  agree  with  the 
prediction  method  of  Volume  II,  Section  3.  The  turboshaft  engine  and  the  TF34 
power  levels  also  agreed  with  this  prediction;  however,  the  low  bypass  ratio 
turbofan  was  6 dB  higher  than  the  levels  predicted  from  high  bypass  ratio 
engine  data.  It  is  suggested  that  the  higher  levels  of  combustor  noise  are 
due  to  lower  turbine  blade  row  attenuation  for  this  low  bypass  engine.  When 
the  blade  row  attenuation  is  accounted  for  by  means  of  a turbine  work  extrac- 
tion term,  the  JT8D-109  data  appear  to  be  within  1 dB  of  the  prediction  line. 

In  addition,  the  Boeing  Company  provided  four  other  sets  of  high  bypass  ratio 
turbofan  engine  data,  which  were  within  ±2  dB  of  the  prediction  line. 

An  evaluation  of  the  Quiet  Engine  Program  Fan  "A"  engine  data  was  per- 
formed over  the  entire  operating  range.  The  data  were  found  to  agree  with 
both  the  overall  power  level  and  the  spectral  shape  used  in  the  prediction 
method  for  power  aettings  up  to  and  Including  takeoff. 

The  combustor  component  spectra  were  found  to  have  the  same  shape  as  the 
prediction  spectrum  and  the  spectral  envelope  (see  Volume  II,  Section  3). 
However,  the  frequency  of  the  peak  level  was  found  to  vary  from  2G0  to  800  Hz; 
whereas  the  prediction  for  engine  core  noise  sets  the  peak  level  at  400  Hz 
± a one-third  octave  band.  It  is  suspected  that  the  influence  of  the  engine 
components  along  the  propagation  path  and  termination  of  the  nozzle  have  an 
Influence  on  both  the  radiated  spectrum  shape  and  peak  frequency,  which  would 
account  for  the  peak  frequency  differences. 

The  overall  power  levels  of  the  component  combustors  were  significantly 
higher  than  the  engine  data  for  any  value  of  the  prediction  parameter.  This 
is  most  probably  due  to  .he  attenuation  occurring  in  passage  through  the 
downstream  turbine  blade  rows.  A theoretical  analysis  for  turbine  blade  row 
attenuation  successfully  models  this  phenomenon,  but  over  predicts  the  atten- 
uation as  indicated  by  the  comparison  of  the  component  and  engine  data. 
Currently,  an  in-depth  program  to  investigate  turbine  blade  row  attenuation  is 
underway  (NAS3-19435  and  DOT-FA75WA-3688) . The  results  of  a parametric  study 
conducted  with  this  analysis  indicated  that  the  low  frequency  noise  attenu- 
ation was  a function  of  the  work  extraction  by  the  turbine  blade  rows.  The 
three  line  prediction  me^  od  for  core  noise  levels  given  in  Volume  II  col- 
lapsed to  a universal  7 ie  when  the  engine  data  were  normalized  by  the  design 
point  work  extraction  by  the  turbines. 
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SECTION  1.0 
INTRODUCTION 


A predict!  ■*.  for  low  frequency  core  engine  noise  was  derived  from  an 
evaluation  of  available  data  as  a part  of  Phase  IV  of  the  Core  Engine  Noise 
Program.  This  basic  (Reference  1)  prediction  was  modified  to  take  into 
account  the  three  haslc  engine  types,  turbofan,  turboshaft  and  turbojet. 

The  correlation  that  was  rhown  was  based  primarily  on  the  low  frequency  data 
from  the  engines  at  low  power  setting. 

The  engine  data  available  at  the  time  was  in  many  cases  sign* ficantly 
influenced  by  jet  noise  at  the  higher  power  settings.  Since  the  writing  of 
that  report,  additional  engine  data  have  become  available  to  enter  into  the 
correlation.  Some  of  these  data  permit  evaluation  of  engine  core  noise  at 
higher  power  settings. 

The  differences  in  the  absolute  level  of  the  low  frequency  noise  radia- 
tion from  the  three  engine  types  has  been  attributed  to  the  varying  attenua- 
tion of  combustor  noise  while  propagating  through  the  downstream  turbine 
blade  rows.  An  analysis  to  explain  the  phenomenon  of  attenuation  of  low 
frequency  noise  by  turbine  blade  rows  has  become  available  and  some  results 
from  that  analysis  are  presented.  The  assumption  of  this  attenuation  mechanism 
permited  addition  of  component  data  to  the  correlation.  Two  series  of 
component  teats  are  reported. 


SECTION  2.0 


ADDITIONAL  ENGINE  DATA 


During  the  evaluation  cf  core  noise  data  presented  in  Volume  III  it  was 
noted  that  several  sets  of  additional  engine  data  were  forthcoming  which  were 
not  available  at  the  time  of  that  report's  publication.  The  engines  were 
the  NASA  JT8D-109  Refan,  a Garrett  turboshaft  and  a TF54  high  bypass  turbofan. 
The  results  of  these  subsequent  engine  tests  are  presented  here.  In  addition, 
the  Boeing  Company  provided  6 data  sets  from  turbofan  engines  for  comparison 
with  the  prediction. 


2.1  NASA  JT8D-109  REFAN  ENGINE 


Low  frequency  core  noise  from  a JT8D-109  refan  engine  has  been  evaluated 
by  comparing  this  data  to  the  core  noise  prediction.*  Low  frequency  core 
noise  is  present  for  this  engine  over  a range  of  lower  power  settings. 

Figure  2-1  shows  that  the  spectral  shapes  fall  between  the  broader  more 
conservative  spectral  envelope  (Volume  III,  Figure  3.3-2)  and  the  more  narrow 
T64  spectrum  which  was  found  to  be  a good  approximation  of  component  confcustor 
noise.  Both  of  these  spectra  are  peaked  at  400  Hz,  (however  the  peak  could  be 
one-third  octave  band  higher  or  lower) . The  Refan  data  peaks  between  400  and 
500  Hz.  The  data  indicates  that  there  is  a variation  in  the  spectrum  shape 
with  angle.  The  refan  spectra  fall  between  the  two  suggested  by  the  prediction 
method,  however. 

The  predicted  directivity  peaks  at  120°,  and  does  not  change  shape  with 
speed.  The  refan  directivities  are  similar  in  shape,  however  the  angle  at 
which  they  peak  changes  from  130°  to  110°  as  engine  speed  increases  (Figure 
2-2) . 


These  variations  in  spectral  shape  and  directivity  are  considered  minor. 
The  major  variation  of  the  JT8D-109  Refan  occurs  in  level.  The  NASA  data  were 
provided  in  the  form  of  peak  0ASFL  versus  percent,  fan  speed,  and  on  this  basis 
the  0ASPL  at  120°  was  found  to  be  6 dB  higher  (Figure  2-3)  than  the  prediction 
for  typical  high  bypass  ratio  turbofan  engines  used  in  the  GE  method.  This 
level  is  almost  as  high  as  the  turboshaft  engines  (which  were  8 dB  higher  than 
tuxbofans) . This  difference  in  levels  could  possibly  be  due  to  either  of  two 
contributing  factors.  The  first  is  that  the  ,,can"-type  combustors  employed  in 
the  JT8D-109  engine  could  possibly  have  gonerated  higher  levels.  The  core 
noise  prediction  was  developed  based  on  annular  combustors.  The  second  factor. 


* These  data  were  provided  by  E.  Krejsa  of  NASA-LEWIS  Research  Center. 
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Figure  2-1.  Conparison  of  Predicted  and  Measured  Spectral  Shapes 
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OASFL  - OASPI 


QNnOS  TTVHSAO 


GE  PREDICTION 
(For  Hi/?  Turbofan) 


70  80  90 

% OF  RATED  FAN  SPEED 


Figure  2-3.  Comparison  of  Predicted  and  Measured  Low  Frequency 

Core  Noise  Levels  for  the  NASA  JT8D-109  Refan  Engine. 
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however,  is  probably  of  greater  significance.  This  Involves  the  work 
extraction  by  the  turbines.  The  transmission  loss  due  to  turbine  blade  row 
attenuation  appears  to  be  a function  of  the  total  temperature  drop  (T4  - T5) 
through  the  turbine  stages  and  the  three  prediction  lines  for  core  noise 
OAFWL  can  be  collapsed  by  normalizing  the  data  through  a 40  log  (T4  - T5) 
factor  - as  is  shown  in  Appendix  A.  It  is  interesting  to  note  that  at  design 
point,  (T4  - T5)  for  the  refan  engine  was  830*  R (461*  K)  while  that  for 
the  typical  high  bypass  turbofan  engine  used  in  the  OAFWL  correlation  was 
1120°  R (622*  K) . The  refan  engine  (T4  - T5)  is  more  representative  of  the 
turboshaft  engines  used  In  the  correlation.  Based  on  these  numbers,  the 
refan  engine  data  would  lie  within  ldB  of  the  correlation  when  normalized 
by  40  log  (T4  - T5)  since  this  would  account  for  5.2  dB  of  the  observed  6 dB 
difference  between  the  refan  data  and  the  high  bypass  engine  correlating 
line. 

2.2  GARRETT  TURBOSHAFT  ENGINE 

The  low  frequency  corn  noise  from  turboshaft  engines  is  relatively  easy 
to  evaluate,  because  fan  noise  does  not  exiBt,  plus  the  jet  noise  is  quite 
low.  A turbo3haft  engine  in  the  840  shaft  horse  power  (SHP)  class  was  tested 
by  Garrett-AlResearch  Company  (Reference  2) . Various  suppression  devices  were 
employed  and  their  effects  measured  with  farfield  microphones  located  every 
15  degrees  on  a 100  ft  (30.48  m)  arc.  The  configuration  evaluated  for  core 
noise  had  a suppressed  inlet  and  an  untreated  eidiaust  duct  9.3  inches  (23.6  cm) 
long.  The  dynamometer  which  absorbed  the  shaft  horsepower  was  sealed  in  an 
acoustical  enclosure. 

A comparison  of  a measured  spectrum  and  the  core  noise  prediction  envelope 
can  be  seen  in  Figure  2-4.  Below  the  first  ground  null,  which  occurs  at  500 
Hz,  there  appears  to  be  the  typical  6 dB  pressure  doubling  due  to  ground 
reflections.  At  the  higher  frequencies,  the  continued  nulls  and  reinforcement 
peaks  which  are  caused  by  the  ground  effect  can  be  seen.  The  levels  of  the 
spectra  which  were  evaluated  are  not  believed  to  be  influenced  by  exhaust  Jet 
noise  since  the  exhaust  velocity  did  not  exceed  250  ft/sec  (76.2  m/sec). 
Therefore,  below  1000  Hz,  the  spectra  are  dominated  by  core  noise.  Turbine 
noise  is  probably  influencing  the  spectra  above  this  frequency.  Keeping  the 
ground  effects  in  mind,  the  prediction  spectrum  peaked  at  400  Hz  is  seen  to 
provide  a good  fit  to  the  data. 

The  radiation  directivity  pattern  for  the  peak  frequency  is  plotted  in 
Figure  2-5.  The  data  scatter  appears  to  be  roughly  what  has  been  seen  from 
the  evaluation  of  other  engine  data.  The  peak  level  appears  to  occur  at  135* 
from  the  inlet  and  is  constant  over  the  speed  range  tested.  The  prediction 
directivity,  which  peaks  at  120°,  still  provides  a reasonable  approximation. 

The  overall  power  levels  were  determined  from  the  turboshaft  engine  power 
level  spectra.  The  prediction  spectra  envelope  level  and  peak  frequency  were 
fitted  to  the  power  level  spectra  as  in  Figure  2-4.  The  overall  level  was 
determined  as  the  peak  level  of  the  spectra  plus  9.9  dB  (where  the  9.9  dB  is 
the  difference  between  the  prediction  spectral  envelope's  peak  level  and  the 
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overall  level  for  that  spectral  shape).  The  overall  power  levels  determined 
by  this  method  are  plotted  against  the  engine  prediction  parameter  in  Figure 
2-6.  The  data  can  be  seen  to  follow  the  same  general  trend  as  the  turboshaft 
engine  prediction  line. 

2.3  NASA  TF34  ENGINE  TEST 

Comparison  of  NASA  data  (Figure  2-7)  for  configurations  with  varying 
amounts  of  treatment  generally  indicate  significant  core  noise  contribution  to 
the  farfleld  levels  for  speeds  below  5100  rpm.  The  low  frequency  noise  suppres- 
sion shown  in  the  spectra  is  due  to  a fairly  short  (12-inch)  and  deep  (4-inch) 
bulk  absorber  used  in  the  core  nozzle  along  with  the  usual  high  frequency 
turbine  treatment. 

The  low  frequency  suppression  cannot  be  seen  at  the  higher  power  settings, 
indicating  jet  noise  domination  of  the  spectrum.  These  observations  were 
verified  by  plotting  the  directivities  for  the  affected  frequencies. 

Earlier  in  the  Core  Engine  Noise  Program,  core  noise  estimates  for  the 
TF34  were  derived  from  General  Electric  data  and  spotted  on  the  OAPWL  predic- 
tion curve.  However,  these  data  Included  considerable  undulations  (possibly 
ground  reflection)  in  the  spectra  and  hence  there  was  some  doubt  whether: 

a.  core  noise  was  present  in  the  farfleld  spectra  (this  question 
has  been  addressed  above),  and 

b.  the  valleys  and  peaks  in  the  spectra  could  be  attributed 
entirely  to  ground  reflections  and,  therefore,  the  freefield 
spectra  defined  by  simply  faring  in  a smooth  line. 

General  Electric  data  are  compared  with  the  NASA  ground  microphone  spectra 
in  Figure  2-8.  A one-to-one  comparison  lo  impossible  since  the  General  Elec- 
tric data  are  for  separate  flow  configurations  (fan  exhaust  upstream  of  core 
exhaust)  while  those  from  NASA  are  for  a confluent  flow  configuration.  It  is 
clear  however,  that  the  undulations  in  the  General  Electric  spectra  (below 
1000  Hz)  do  result  from  ground  reflections.  The  agreement  between  the  General 
Electric  and  NASA  levels  indicate  that  the  early  General  Electric  estimates  of 
TF34  core  noise  were  fairly  good. 

2.4  BOEING  DATA 


The  Boeing  Company  has  provided  data  from  six  turbofan  engines  in  a 
comparison  with  the  prediction.  The  overall  sound  pressure  levels  at  the  120° 
angle  on  a 200  foot  (60.96m)  sideline  were  plotted  (Figure  2-9)  against  the 
prediction  parameter.  It  was  observed  that  the  high  bypass  ratio  turbofan 
engine  data  agreed  very  well  with  the  turbofan  prediction  line.  However,  the 
low  bypass  turbofan  fell  close  to  the  turboshaft  1j  ie  similar  to  the  results 
(Figure  2-3)  for  the  Refan. 

Note  that  the  data  in  Figure  2-9  extend  to  power  settings  beyond  the 
original  correlation. 
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Figure  2-7.  TF34  Spectral  Comparison  Using  Engine  Center-Line  Microphones 


Figure  2-9.  Comparison  of  Boeing  Data  with  GE  Procedure 


SECTION  3.0 


ENGINE  HIGH  POWER  SETTING  EVALUATION 


Examination  of  acoustic  data  from  turbojets  and  turbofans  at  high  power 
settings  revealed  that  in  most  cases  jet  noise  contaminated  or  overshadowed  the 
core  noise  in  the  aft  quadrant.  In  order  to  provide  additional  high  power 
setting  data  an  Investigative  technique  consisting  of  studying  the  exponent 
of  the  exhaust  velocity  for  each  frequency  at  each  angle  along  with  the  effects 
of  acoustic  treatment  on  the  contaminating  noise  sources  was  used. 

Acoustic  data  was  evaluated  for  discernible  core  noise  at  hxgh  power 
settings  on  the  T64,  GE12,  Quiet  Engines  "A"  and  "C",  TF34,  CF6 , J65,  J79, 

J85  and  GE4. 

There  was  no  problem  with  jet  or  fan  noise  contamination  for  the  turbo- 
shaft  engines  (T64,  GE12) . Data  encompassing  the  entire  operating  range  is 
spotted  on  the  prediction  line.  The  Garrett  data  also  included  high  power 
settings . The  examination  of  turbofan  engine  data  revealed  that  Engine  "A" 
data  provided  discernible  core  noise  levels  over  the  entire  operating  range  of 
the  engine.  The  data  from  the  other  engines  mentioned  above  did  not  yield  any 
significant  or  consistent  indications  of  core  noise  at  high  power  settings. 


3.1  THE  SPECTRUM  LEVEL 

The  primary  problem  of  evaluating  core  noise  at  the  higher  power  settings 
is  the  contaminating  influence  of  other  noise  sources,  mainly  the  fan  and 
the  jet.  There  are  two  methods  that  can  be  used  to  evaluate  the  relative 
effects  of  these  other  noise  sourcee.  The  effect  of  fan  noise  can  be  evalu- 
ated by  comparing  the  spectral  data  from  an  unsuppressed  fan  to  a fully 
suppressed  fan  configuration  at  the  angle  where  core  noise  is  expected  to 
peak,  namely  120°  from  the  engine  inlet  (see  Section  3,  Volume  III).  When 
the  levels  of  the  two  spectrums  converge  it  can  be  reasonably  assumed  that  the 
fan  has  little  or  no  effect  on  the  spectrum. 

The  presence  of  jet  noise  can  be  determined  by  examining  the  sound  pres- 
sure levels  for  the  individual  one-third  octave  frequency  bands  and  plotting 
them  against  the  core  exhaust  velocity  over  the  operating  range  of  the  engine. 
When  the  slope  of  the  plot  indicates  an  exhaust  velocity  to  the  sixth  power 
or  higher  relationship,  it  can  be  assumed  that  jet  noise  is  dominating  the 
spectrum.  However,  when  the  slope  of  the  data  indicates  an  exhaust  velocity 
to  the  fourth  power  dependency,  it  can  be  assumed  that  some  other  source  domi- 
nates. Hence,  the  one-third  octave  band  levels,  which  could  not  be  identified 
as  either  fan  noise  or  jet  noise  and  which  followed  a core  exhaust  velocity 
to  the  fourth  power,  are  indicative  of  core  noise. 
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Data  from  Engine  "A"  of  the  NASA  Quiet  Engine  Program  provide  a good  set 
of  data  for  evaluation  of  core  noise  at  the  higher  power  setting.  Figure  3-1 
shows  a low  power  setting  spectrum  at  120°  from  the  inlet.  The  region  of 
fan  noise  dominance  can  be  seen  beginning  at  1250  Hz  and  continues  to  10,000 
Uz.  On  the  very  low  frequency  end,  jet  noise  can  be  seen  to  dominate.  The 
prediction  for  coannular  jet  noise  is  plotted  over  the  measured  data.  The 
remaining  portion  of  the  spectrum,  between  400  and  1000  Hz,  has  the  core  noise 
prediction  spectrum  drawn  in.  The  influence  of  jet  noise  at  any  particular 
one  third  octave  band  is  evaluated  from  the  exponent  of  the  exhaust  velocity 
(Figure  3-2)  . At  315  Hz  a line  representing  the  exhaust  velocity  raised  to 
the  sixth  power  has  been  drawn  through  the  data.  The  close  fit  of  the  data  at 
this  frequency  indicates  that  jet  noise  is  dominating  over  the  whole  engine 
speed  range  with  the  exception  of  the  lowest  speed  point.  This  lowest  point 
diverges  in  a direction  which  would  indicate  a velocity  to  the  fourth  power, 
or  the  presence  of  core  noise.  This  same  trend  is  observed  at  400  Hz.  At 
500  Hz  however,  there  is  a divergence  from  the  velocity  to  the  sixth  power 
which  now  occurs  at  a higher  power  setting.  At  this  frequency  the  lowest 
three  points  appear  to  be  core  noise  and  the  upper  3 points  appear  to  be  jet 
noise  dominated.  The  next  three  frequencies,  (630,  800  and  1000  Hz)  clearly 
follow  the  exhaust  velocity  to  the  fourth  power  relationship  over  the  entire 
engine  speed  range. 

The  frequency  range  of  core  noise  dominance  as  a function  of  engine 
power-level  was  evaluated  from  the  spectra.  Figure  3-3  shows  three  engine 
spectra,  each  with  the  component  sources  which  generate  the  spectrum  indicated. 
The  level  of  the  core  noise  can  be  determined  by  fitting  the  core  noise  pre- 
diction spectrum  to  the  engine  data.  Based  on  the  peak  level  of  this  spectrum 
and  the  core  noise  prediction  directivity  the  core  noise  power  levels  can  be 
obtained  over  the  engine's  entire  speed  range.  The  power  levels  calculate  by 
this  method  are  plotted  against  the  core  noise  prediction  parameter  in  Figure 
3-4.  As  can  be  seen,  the  level  follows  the  prediction  line  for  turbofan 
engines  very  well. 

These  data  verify  that  the  core  noise  prediction  parameter  is  valid  over 
the  whole  operating  range  of  an  engine. 


3.2  SPECTRAL  SHAPE  AT  HIGH  POWER  SETTINGS 

At  the  lower  power  settings,  the  core  noise  prediction  spectrum  fits  the 
measured  Engine  "A"  data  for  6 one-third  octave  bands  (Figure  3-3) , When 
the  engine  epeed  is  increased  to  95%  fan  speed  the  slope  of  the  jet  noise 
spectrum  approaches  the  slope  of  the  core  noise  spectrum;  however  the  core 
noise  spectrum  still  maintains  a good  fit  over  5 one-third  octave  bands.  There 
is  no  indication  that  the  core:  noise  spectral  shape  differs  appreciably  from 
the  prediction  spectrum  over  the  engine  operating  range. 
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3.3  ENGINE  HIGH  POWER  SETTING  DIRECTIVITY 


At  the  higher  power  settings  It  was  found  that  for  angles  greater  than 
120°  the  directivity  pattern  followed  the  jet  noise  prediction.  At  .an  speeds 
higher  than  approximately  70X , jet  noise  has  a significant  influence  on  the 
directivity.  In  this  case  it  is  assumed  that  the  core  nolle  directivity  does 
not  change  at  the  higher  pcwer  levels. 

As  was  noted  earlier,  the  turboshaft  data  were  found  to  include  cignifi- 
cant  core  noise  over  the  entire  operating  range,  as  were  those  from  Engine 
"A".  The  data  from  the  other  engines  were  found  to  be  controlled  by  either 
jet  or  fan  noise. 

The  Boeing  data  (Figure  2-9)  also  extends  to  high  power  settings  and 
thus  provide  additional  verification  of  the  prediction  line. 
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SECTION  4.0 


ADDITIONAL  COMPONENT  DATA 


4.1  COMPONENT  COMBUSTOR  SPECTRA 


Component  combustor  tests  were  preformed  at  atmospheric  pressure  under 
the  original  portion  of  this  program.  The  results  are  reported  in  Volume  II 
(Section  3) . The  acoustic  data  were  originally  acquired  with  16  ft  (4.88  m) 
high  microphones  on  the  40  ft  (12.19  m)  radius  arc,  however,  ground  reflections 
were  found  to  create  nulls  and  peaks  in  the  low  frequency  rogion  of  the  spectra, 
A spectral  correction  for  these  ground  effects  was  obtained  from  a comparative 
examination  of  data  taken  with  microphones  at  the  16  ft  (4.88  m)  and  2 ft 
(0.61  m)  height:.  The  two  foot  microphone  height  moves  the  first  null  of  the 
ground  reflection  to  a frequency  of  about  2000  Hz  and  out  of  the  range  of 
interest.  These  corrections  were  applied  to  the  original  data  taken  at  the 
higher  microphone  height.  The  resulting  spectra  are  plotted  in  Figures  4-1 
and  4-2  for  the  Advanced  Technology  (AT)  and  CF6  combustors  respectively.  For 
comparison  purposes  the  core  noise  prediction  spectrum  shape  is  shown  for  one 
temperature  setting.  The  prediction  spectrum  shape  agrees  with  the  shape  of 
the  data,  however  the  frequency  at  which  the  peak  occurs  is  200  Hz.  This  is 
two  one-third  octave  bands  lower  'han  what  would  be  expected  from  the  current 
engine  core  noise  prediction  procedure.  The  difference  may  be  due  to  the 
downstream  turbine  blade  rows  and  the  nozzle  end  impedance  encountered  in 
engines.  Cross  sectional  schematics  of  a turbofan  engine,  an  atmospheric 
pressure  component  combustor  test  rig,  and  a high  pressure  component  combustor 
test  configuration  are  shown  in  Figure  4-3. 

Noise  measurements  of  component  combustors  operating  at  high  inlet  tem- 
peratures and  pressures  were  taken  by  GE  under  the  Experimental  Clean  Combustor 
Program  (Contract  No.  NAS3-18551) . This  program  involved  testing  three 
design  philosophies  of  low  emission  full  scale  annular  conbustors,  the  Swirl- 
Can,  Radial /Axial,  and  Double  Annular.  The  combustors  were  tested  in  a high 
pressure  component  test  rig.  Acoustic  probes  measured  the  levels  immediately 
downstream  of  the  combustor.  Representative  sound  pressure  level  spectra  can 
be  seen  in  Figure  4-4.  Simulated  approach  and  takeoff  combustor  inlet 
conditions  are  shown.  The  spectral  envelope  peaks  at  500  Hz  and  is  seen  to  be 
in  good  agreement  with  the  measured  spectra  for  all  the  configurations  at 
approach  power  conditions.  The  T64  spectrum,  which  falls  off  somewhat  rapidly 
above  1000  Hz,  is  also  shown.  At  the  takeoff  (inlet)  condition,  however, 
some  of  the  configuration's  peak  levels  appear  to  shift  to  a higher  frequency. 
The  spectral  envelope  is  peaked  at  800  Hz  and  fits  the  Double  Annular  and  the 
90-Swirl-Can  spectra.  The  Radial/Axial  combustor  shows  higher  levels  at  the 
higher  frequencies,  but  can  not  be  considered  a typical  annular  combustor. 

This  configuration  employed  afterburner  type  f lameholders  which  are  believed 
to  be  responsible  for  the  broader  spectral  shape.  Again,  the  end  termination 
was  somewhat  different  from  that  for  engines,  and  this  may  involve  a spectral 
re-ordering  due  to  frequency  dependent  attenuation  and  radiation. 
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FIGURE  4-1  POWER  LEVEL  (CORRECTED  FROM  16  FOOT  HIGH  TO 

2 FOOT  HIGH  MICRpPHQNES) , F101  COMPONENT  COMBUSTOR 
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Figure  4-2.  Power  Level  (Corrected  from  16  foot  High  to 

2 foot  High  Microphones)  C76  Component  Combustor. 
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The  higher  pressure  combustor  noise  tests  appear  to  fit  the  broader 
spectral  envelope,  while  the  engine  data  in  general  tend  to  follow  the  T6A 
(narrower)  spectrum.  The  difference  may  be  attributed  to  the  effect  of  the 
turbine  blade  rows  downstream  of  the  combustor  in  the  case  of  engines. 


A. 2 COMPONENT  COMBUSTOR  POWER  LEVEL 


The  noise  levels  of  the  component  combustors  can  be  compared  on  a power 
level  basis.  The  CF6  and  AT  combustor  power  levels  were  calculated  from  far- 
field  arc  data.  Their  overall  levels  were  arrived  at  by  logarithmically 
summing  the  one-third  octave  band  power  levels  over  the  frequency  range  of 
100  to  2000  Hz.  The  high  pressure  combustor  data  were  summed  over  the  31.5  to 
5000  Hz  frequency  range  and  the  power  level  was  calculated  based  on  the  local 
temperature,  pressure,  Mach  number  and  cross  section  annulus  area: 


PWL 


SPL  + 10Log1Q  (A)  + 10Log10 


+ 10Log10 


where  A is  the  annulus  area  (square  meters)  at  the  probe's  sensing  holes,  P 
and  T are  the  local  pressure  (atmospheres)  and  temperature  (°K)  where  the 
subscripts  indicate  standard  day  reference  levels  and  M^  is  the  Mach  number  of 
the  flow  past  the  probe.  The  constant  K is  the  sum  of  the  SPL  and  PWL  refer- 
ence levels  and  the  standard  day  characteristic  impedance  of  the  medium.  It 
should  be  noted  that  there  is  no  coefficient  in  front  of  the  Mach  number.  A 
coefficient  may  be  used  if  it  is  assumed  the  acoustic  waves  strike  the  probe 
at  random  incidence  angles.  The  wavelength  of  the  sound  being  considered  is 
much  larger  than  the  passage  height  dimension,  therefore  a plane  wave  assump- 
tion can  be  made  which  makes  the  coefficient  unnecessary.  The  standard  day 
reference  temperature  and  pressure  are  288°  K and  one  atmosphere.  The  constant 
K is  9.9  dB. 

These  overall  power  levels  are  compared  by  plotting  them  against  the 
engine  prediction  parameter  (Figure  A-5)  . The  major  portion  of  the  component 
data  fall  above  the  turbojet  prediction  line. 

The  grouping  of  the  3 classes  of  engines  and  the  relationship  of  the 
component  noise  power  levels  to  the  engine  data  strongly  suggests  an  attenua- 
tion mechanism  for  combustor  noise  along  the  propagation  path  through  the 
turbine  stages  and  exhaust  duct. 


A. 3 TURBINE  BLADE  ROW  ATTENUATION 


The  attenuation  of  low  frequency  noise  through  a turbine  stage  has  been 
analytically  predicted  by  K.  Bekofske  (Reference  3).  His  study  is  based  on 
an  actuator  disk  model.  An  actuator  disk  is  a cascade  model  with  infinitesi- 
mal chord  length  and  blade  spacing.  Since  combustor  low  frequency  acoustic 
energy  in  the  high  temperature  environment  has  characteristic  wavelengths  that 
are  large  when  compared  to  both  the  chord  and  transverse  spacing  of  a blade 
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Figure  4-5.  Component  Combustor  Power  Level  Vs.  Prediction  Parameter. 


row,  the  use  of  the  actuator  disk  model  Is  justified.  In  this  analysis,  the 
relative  Mach  nuaber  at  both  the  inlet  and  discharge  side  of  either  the  nozsle 
or  bucket  has  been  assumed  to  be  subsonic. 

lhe  stage  attenuation  is  the  difference  between  the  initial  wave  and  the 
wave  amplitude  transmitted  through  the  stage.  This  difference  is  determined 
by  a reflection  coefficient.  An  infinite  duct  is  assumed  upstream  of  the  blade 
row  such  that  the  reflected  ray  continues  propagating  upstream  and  is  completely 
removed  from  the  system.  Hie  calculation  of  this  reflection  coefficient  takes 
into  effect  the  changes  in  Mach  number  and  the  flow  turning  angle.  The  stage 
attenuation  can  be  calculated  as  a function  of  the  initial  wave  incidence 
angle.  It  is  assumed  that  the  initial  wave  incidence  angle  on  each  stage  is 
close  to  zero  degrees,  (which  represents  a plane  wave  striking  the  stage 
normally.  Figure  4-6)  . In  actuality  it  is  expected  that  there  is  some  distri- 
bution of  the  incidence  wave  as  a function  of  angle,  with  the  highest  prob- 
ability at  zero  degrees. 

The  input  to  a computer  program  based  on  this  analysis  (listed  in  Appendix 
B)  consists  of  the  axial  velocity,  absolute  flow  angle,  static  pressure  and 
total  temperature  at  three  locations  in  any  given  stage:  upstream  of  the 

nozzle,  between  the  nozzle  and  bucket,  and  downstream  of  the  bucket.  The 
aerodynamic  parameters  of  each  stage  of  both  the  high  and  low  pressure  turbines 
of  a typical  high  bypass  turbofan  engine  for  a takeoff  condition  were  input  to 
this  computer  program.  Figures  4-7  and  4-8  show  the  first  stage  of  the  high 
pressure  turbine  and  the  last  stage  of  the  low  pressure  turbine  for  a typical 
high  bypass  ratio  engine.  Both  figures  show  an  angle  where  the  attenuation 
sharply  Increases.  Incident  waves  whose  angles  are  larger  than  these  angles 
will  be  cut  off  and  theoretically  no  portion  of  their  acoustic  energy  will  be 
transmitted. 

Table  4-1  lists  the  stage  total  attenuation  at  zero  degrees  for  each 
stage  and  the  total  attenuation  for  the  seven  turbine  stages.  A total  of 
38.9  dB  is  predicted  for  this  takeoff  condition.  From  the  examination  of  the 
low  pressure  turbine  (Table  4-2)  attenuation  as  a function  of  engine  speed,  it 
was  found  that  the  attenuation  does  not  appreciably  change  with  power  setting. 

A comparison  of  the  component  and  turbofan  correlating  line  Indicates  that  the 
analysis  somewhat  over  predicts  the  attenuation  occurlng  in  transmission  through 
the  blade  rows.  This  over  prediction  of  the  attenuation  might  be  attributed 
to:  (1)  the  assumption  that  the  energy  reflected  at  a blade  row  is  completely 

removed  from  the  system.  In  reality,  the  reflected  wave  can  propagate  upstream 
to  the  next  blade  row  and  part  of  this  may  be  reflected  back  downstream, 
effectively  reducing  the  attenuation  and  (2)  conditions  in  an  engine  which 
deviate  from  the  assumed  mean  flow  between  stages.  In  particular,  the  inci- 
dence angle  may  vary  from  zero.  The  trend  in  the  measured  data  and  the 
analytic  prediction  nevertheless,  both  indicate  significant  attenuation  of 
noise  propagating  through  the  turbine.  Further  development  of  the  blade  row 
attenuation  prediction  technique  Is  required.  This  work  has  been  Initiated 
under  Contract  NAS3- 19435. 
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Table  4-1.  Low  Frequency  Nolee  Attenuation  for 
High  Bypass  Turbofan  Engine  Turbines 


Attenuation 

Stage 

dB  at  cero  degraea 

1 

11.3 

High 

2 

9.2 

Pressure 

Turbine 

1 

5.2 

2 

4.2 

3 

4.0 

Low 

4 

3.0 

Pressure 

5 

2.0 

Turbine 

38.9  dB  Total 

Table  4-2. 

Variation  of  Low  Frequency  Noise 

Attenuation  with  Speed. 

Low  Pressure  Turbine 

Attenuation 

X Fan  RPM 

Stage 

dB  at  aero  degrees 

1 

5.0 

2 

3.8 

91 

3 

3.75 

4 

2.8 

5 

1.8 

17.15 

1 

5.25 

2 

4.1 

88 

3 

3.8 

4 

2.8 

5 

1.9 

17.85 

1 

5.0 

2 

3.8 

^6 

3 

3 . 5 Approach 

4 

3.5  Power 

5 

1.5 

17.3 

1 

4.9 

2 

3.6 

65 

3 

3.2 

4 

2.2 

5 

1.3 

15.2 
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SECTION  5.0 


CONCLUSIONS 


It  is  evident  that  more  needs  to  be  understood  about  the  nature  of 
the  combustor  noise  source  and  the  propagation  path  it  follows  through 
an  engine  before  this  prediction  procedure  can  be  refined  further.  Never 
theless,  this  prediction  technique  is  quite  good  when  the  range  of  air 
flow  rates  (1  to  600  lb/sec)  of  the  engine  data  which  have  been  evaluated 
is  considered.  The  important  points  in  this  report  are  tabulated  here: 

e The  difference  between  component  combustor  levels  and  engine 
levels  indicates  attenuation  along  the  propagation  path  through 
the  engine. 

e The  overall  power  level  prediction  procedure  using  three  lines 
for  the  engine  classes  is  valid  for  the  entire  operating  range. 
The  data  can  be  collapsed  into  a single  line  correlation  by  in- 
cluding a turbine  work  extraction  term  to  account  for  the  blade 
row  attenuation. 

e The  spectral  shapes  of  combustor  and  engine  core  noise  vary  some- 
what but  fall  between  the  prediction  spectrum  and  the  spectral 
envelope. 

e The  frequency  at  which  the  spectrum  peaks  also  varies,  however, 
based  on  engine  data  the  peak  level  occurs  at  400  Hz  + a one- 
third  octave  band. 

e The  radiation  directivity  pattern  was  generally  found  to  peak 
at  120°  from  the  engine  inlet. 
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APPENDIX  A 


THE  ENGINE  LOW  FREQUENCY  CORE  NOISE  PREDICTION 


The  procedure  consists  of  first  calculating  the  overall  power  level 
(OAPWL)  from  the  prediction  parameter.  Proa  that,  the  one-third  octave  band 
power  levels  are  determined  by  iaposing  the  prediction  spectrum.  These  power 
levels,  the  arc  radius,  and  the  prediction  directivity  are  finally  used  to 
calculate  the  one-third  octave  band  sound  pressure  levels  for  any  angle  around 
the  arc. 

The  overall  power  level  correlation  based  on  engine  data  was  observed  to 
form  three  separate  but  parallel  prediction  lines  representing  the  three 
engine  types:  turbojet,  turboshaft,  and  turbofan.  The  analysis  for  the 

attenuation  of  low  frequency  by  turbine  stages  suggests  that  the  work  extrac- 
tion by  the  turbines  is  a significant  parameter(4)  in  this  delineation  of  the 
engine  types.  The  temperature  drop  across  the  turbines  at  design  point  was 
introduced  into  the  correlation  and  a data  collapse  achieved  as  shown  in 
Figure  A-l.  The  prediction  line  is  defined  by: 

OAPWL  - 171  + 10  log10  [*  (T4-T3)2(p3/p0)2]  - 10  log  (TA-T5)pe8  (A-l) 


where  W is  the  weight  flow  rate  in  lb/sec. , T denotes  the  total  temperature  in 
degrees  Rankine,  and  p the  density;  station  0 is  ambient,  3 the  combustor 
inlet,  4 the  combustor  exit,  and  5 the  low  pressure  turbine  exit.  The  design 
point  (T4-T5)  must  be  used  for  any  engine,  regardless  of  the  power  setting. 

Actually,  the  best  mean  square  fit  to  the  data  of  Figure  A-l  is  provided 
by  a line  with  slope  of  1.05  instead  of  the  slope  of  1.00  as  used  in  equation 
(A-l).  However,  the  difference  between  the  two  lines  is  very  small,  about  1 
dB  at  either  extreme,  and  the  standard  deviation  is  about  2.6  for  both.  Since 
either  line  is  well  within  the  scatter  range  of  the  other,  the  final  choice 
was  made  on  the  basis  of  understanding  the  causation:  the  prediction  param- 

eter being  a product  of  the  heat  release  rate  [10  log  ti(T4-T3)]  and  the 
thermo-acoustic  efficiency  [10  log  (T4-T3)]. 

The  temperature  drop  term  is  related  to  the  turbine  pressure  ratio  which 
has  previously  been  used  in  correlations  by  Gerend,  et  al.(5)  and  Grande^). 
The  number  of  turbine  stages  cannot  be  recommended  for  use  in  a correlation  as 
the  analysis  strongly  indicates  that  this  is  not  a linear  dependency  (see 
Table  4-1).  The  effects  of  the  off-design  performance  on  the  turbine  attenu- 
ation and  nozzle  termination  are  as  yet  largely  unknown  factors.  Some  inves- 
tigation of  thes  effects  is  being  conducted  under  contracts  NAS3-19435  and 
DOT-FA75WA-3688.  Meanwhile,  equation  (A-l)  is  recommended  for  predicting  the 
OAPWL. 
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The  peak  level  of  the  prediction  spectrum  (Figure  A-2)  is  given 
by: 

Peak  PWL  = OAPWL  - 0.8,  dB  (A-2) 

This  sets  the  level  of  the  peak  frequency  one-third  octave  band  power 
spectrum.  The  peak  frequency  is  nominally  set  at  400Hz,  although  it 
is  realized  that  it  may  occur  at  a one-third  octave  band  to  either  side, 
that  is  at  315  or  500  Hz. 

The  complete  power  level  spectrum  is  defined  by  utilizing  the  narrow 
(T64)  prediction  spectrum  of  Figure  A-2  with  the  peak  PVL.  However,  if 
a degree  of  conservatism  is  desired,  the  broader  spectrum  envelope  shown 
in  Figure  A-2  can  be  used  insteady. 

The  sound  pressure  level  at  any  point  around  the  arc  can  be  found 
from  the  power  levels  through: 

SPL  = PWL  - 20  Log  R - 9.25  + DI  - AA  (A-3) 

Where  R is  the  arc  radius  in  feet  and  DI  is  the  directivity  index  from 
Figure  A-3,  and  AA  is  the  air  attenuation. 
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APPENDIX  B 


PROGRAM  LISTING  - PREDICTION  OF  BLADE  ROW  ATTENUATION 


ELI ST  TAD** 
12/11/73  11137 


00010  REAL  Kl#  K2»  KRA#  KX0VKY#  Ml#  M2#  MIA#  M2A#  Ml T#  M2T 
00020  REAL  M3A#  M4A#  M3T#  M4T 

00030  DIMENSION  Q9C1005#  ATTEM2C 1 00 > » QOSAVEOOO) 

00040  P I “3 • 1 41 5927  i GRA V=32 • 1 7 40 5 1C0NV*778.26 
00050  GAMMA* 1 • 32  1RBC*0. 06854 

00060  PR I NTlV ABSOLUTE  VELOCITIES  ARE  INPUT  IN  OPTION  1" 

00070  PRI NTt "R ELAT I VE  MACH  NUMBERS  ARF.  INPUT  IN  OPTION  2" 

00080  PRI  NT!  "SELECT  OPTION  1 0R  2 " IREADt  INPUT 

00090  I FC INPUT  .EO.  2)  G0  T0  130 

00100  PRINTl "ABSOLUTE  VELOCITIES  ARE  IN  FEET/SECOND"  1PRINT  1460 
00110  PRINTl "ABSOLUTE  FL0W  ANGLES  ARE  IN  DEGREES  RELATIVE  T0  PACE  PL  ,INE" 
00120  PRINTl"  AND  NEGATIVE  IF  FLOW  DIRECTION  OPPOSITE  WHEEL  ROTATION" 
00130  130  PRINT  1460  I PR  I NTl  "PR  ESSURES  ARE  IN  PSl'A" 

001  40  PRINT  1460  1 PR I NTl "TEMP ERATURES  ARE  IN  DEGREES  RANKINE" 

00150  150  N*0  J PR  I NT  1460 
00150  DO  170  I » 1 > 1 00#  1 
00170  170  ATTEN2CI>=»0. 

00180  I FC INPUT  • E9  • 25  GO  TO  290 

00190  P R I N T 1 " I NP UT  DATA  UPSTREAM  OF  NOZZLE" 

00200  PRINTl  "AXIAL  VELOCITY#  ABSOLUTE  FLOW  ANGLE"  I READ!  VIA#PSI^ 

00210  PRI NTl "STATIC  PRESSURE#  TOTAL  TEMPERATURE"  IREADt  PSl* TZ1 
00220  PRI NTl " I NPUT  DATA  BETWEEN  NOZZLE  AND  BUCKET" 

00230  PRI NTl "AXIAL  VELOCITY#  ABSOLUTE  FLOW  ANGLE"  1 READ!  V2A#PSI2 
00240  PRI  NTl  "STATIC  PRESSURE#  TOTAL  TEMPERATURE"  IREAlii  PS2#  TZ2 
00250  PR  I NTl  "INPUT  DATA  DOWNSTREAM  OF  BUCKET" 

00260  PRINT, ("AXIAL  VELOCITY#  ABSOLUTE  FLOW  ANGLE"  IREADt  V3A#  PSI3 
00270  PRI NTt "STATIC  PRESSURE#  TOTAL  TEMPERATURE"  1 READ! PS3*  TZ3 
00230  PRINTl  "INPUT  WHEEL  SPEEDCFPS!"  IREAD!  UW23  !UW*0. 

00290  290  I F< INPUT  .EQ.  1)  GO  TO  420 
00300  PRI NTl "INPUT  DATA  UPSTREAM  OF  N0ZZLE" 

00310  PRI NTt"AXIAL  AND  TANGENTIAL  RELATIVE  MACH  NUMBERS"! READ! Ml A# Ml T 
00320  PRINTl"STATIC  TEMPERATURE  AND  PRESSURE"  1 READ!  TEMP!#PS1 
00330  PRINT! "INPUT  DATA  DOWNSTREAM  OF  NOZZLE" 

00340  PRINTl"AXlAL  AND  TANGENTIAL  RELATIVE  MACH  NUMBERS"! READ! M2A# M2 T 
00350  PRINT! "STATIC  TEMPERATURE  AND  PRESSURE"  IREADt  TEMP2»PS2 
00360  PRI  NTt  " I NPUT  DATA  UPSTREAM  OF  BUCKET" 

00370  PRINTf’AXIAL  AND  ^TANGENTIAL  RELATIVE  MACH  NUMBERS"! READ! M3A# M3 T 
00330  TEMP3*TEMP2  JPS3=PS2 

00390  PRINTflNPlJT  DATA  DOWNSTREAM  OF  BUCKET" 

00,400  PR  I N Tl  "AX  I AL  AND  TANGENTIAL  RELATIVE  MACH  NUMBERS"!  READ!  M4A#  M4T 
00410  PRINT!"STATIC  TEMPERATURE  AND  PRESSURE"  iREADt  TEMP4#PS4 
00420  420  IF< INPUT  . E9.  25  GO  TO  440 

00433  PSI1*PI*PSII/180.  IPSI2*PI*PSI2/180.  ! PS I 3=PI *PSI 3/1 80 . 

00440  443  N=N+1  I I FC  INPUT  • E3«  2)  GO  TO  530 
00450  V1T*V1A/TANCPSI»  > 1 V2 T»V2A/ TANCPSI2) 

00460  V1=SQRTCV1A*V1A+-VI  T*V1T,5  ! V2  = S0RTC  V2A*  V2A*  V2T*  V2T5 

00470  W 1 A = V 1 A JW>A  = V2A  !WIT”V1T-UW  !W2r=V2T-UW 

00430  Wl  = S9RT<  W1  A + WIA+Wl  T*W)T>  1 W2=>S9RTC  W2A*W2A«-W2T*W2T> 

00490  AZ 1 * SORT  l GAMMA+GRA  V*RG(!*CONV*TZ  1 ) 

00503  AZ2*S0RT<  GAMMA  KGRAVt>RGCtC0NV*TZ2) 
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onst  r»  TEMPI «TZl*(  1 .-(GAMMA- l . 3*<  VI  A*  VI  A*  VI  T*Vl  T> /(  2 . *AZ  I *AZ1  ) ) 
oos?3  TFM»2aTZ2*(  i .-(gamma-i  .i*<  v2a*v2a*  V2T+ V2  n/(2.*AZ2+AZ2>i 
00  )30  500  Al=SORr(GAMMA*GRAV*RGC*C0NV*TEMPl  ) 

00 5 <10  A:»  = S9R  r(GAMM4*GRAV*RGC*C0NV*  TEMP33 
00550  ARA=AI/A2  »KR4=1./ARA 
00560  I r < I vjptj  r • EQ • 2)  G0  T0  580 

00570  MIA=W14/AI  3 M2A=W2A/A2  3MIT3W1T/A1  3M2T3W2T/A2 

90  550  550  'll  aS3RT(M1  A*M1  A + Ml  T*M1  T)  I M2=SQR T( M2A*M2A+M2T*M2T> 

00590  0TCt=ArAN2(SQRT(l.-M2A*M2A3,  -M2A) 

09600  9T023A rAN2( -SQRT( 1 • -M2A*M2A3,  -M2A3 

00  610  GTO  I =S  1 'KOTO  1 3/(KRA+<  1 . ♦ M2A*C0S(  OTC  1 » +M2TKSI N<  0 TC 1 ) > ) 

00  620  GT02  = S 1 N(  0 TC2  ) /(  KRA*(  » • +M2A*C0S( Q TC2  3 *M2T*SI  N(  QTC2)  ) » 

00630  GTC3GTC 1 3QTC=QTCI  3M*0 

00640  640  M3M+1  3 QQTC 3 1 80 . *Q TC /P I 

00650  A r03(  1 . -GTC*M1  T3*(  1 .-GTC*Mt  T>  -GTC*G  rC 

00560  OTC3-2.*GTO*M1 A*< 1 . -GTC*M1 T3 

00670  C TC  = - G rC*GTC*(  1 • -.'ll  A*M1A3 

00680  X = 8TC*BTC“4.*A  TC*CTC 

00690  I F( 6 .LT.  0.)  G0  T0  820 

00700  ODTC3( -BTC-SIGN( l .,BTC3+SQRT(BTC*BTC-4.*ATC*CTC33/2. 

00710  QC3ATAN2(DDTC,ATCI  3 QQC3 1 80 . *QC/P I 3PRINT  1460 
00720  IF(N  • EQ  • 2)  G0  T0  740 

007:30  PRI  NT:  “NOZZLE  ATTENUATION**  3PR1NT  1460  3 G0  T0  750 

00740  740  PRI NTS  "BUCKET  AND  STAGE  ATTENUATION"  3 PRINT  1460 

00750  750  PR  I NT:  "UPSTREAM  CUT-0FF  OCCURS  AT  INCIDENCE  ANGLES  0F" 

00760  CUT0FF=9O.+18O.*ARSINCM1A3/PI  3 PRINT  1430/  CUTOFF 

00770  CUT0FF=-CUT0FF  3 PR I NT  1430/  CUT0FF  3PRINT  1460 

00780  PRINT: "DOWNSTREAM  CUT-0FF  OCCURS  AT  A TRANSMISSION  ANGLE  OF" 

00790  PRINT  1430,  QQTC  3 PRINT  1460 

00800  PRINT: "THE  CORRESPONDI NG  INCIDENCE  ANGLE  IS" 

00810  PRI  NT  1 430,  QQC  3 PRINT  1460 
00820  520  I F( M .EQ.  23  GO  TO  840  ' 

00830  GTC  = GTC2  t 9TC3QTC2  3 G0  T0  6*40 
00840  840  1 F( Ml  .LT.  1.  .AND.  M2  .LT.  1.3  G0  T0  860 
00550  PRI NT*  "SU°ERS0NIC  RELATI  VE  s-LOW"  3G0  T0  1400 
00860  860  I F( N . EQ . I 3 GO  T0  890 

00870  PRINT:"  THETA  THETA  T T V R 

00830  A ATTEN  ATTEN2  THETA  I"  3.  GO  TO  900 

00390  890  PRINTS"  THETA  I THETA  T T V R ATTEN” 

09900  900  RH0RAaP31 *TEMP2/(PS2*TEMP1 3 

00910  ALPHA  = ATAN2<  Ml T/Ml A3  3 BETA* ATANSC M2T, M2A 3 

099:29  G0  T0  (930,  19593,  N 

00930  930  99(  1 3 a90  . ♦ 1 30  . *ARS  I N<  Ml  A 3 /PI  3 QQMI  N3 -QJ9C  I 3 3S310. 

00949  9QH0LD=1O.*4INT( CQQ(1 3 + 10. 3/10.3 
09950  I F( QQC  .LT.  0.3  30  T0  980 

00960  QQ(  1 3=-90 . -1 80 • +ARSI N<  Ml  A3 /P I 3 QQMAX3?QQ( 1 3 3S3-10. 

09970  QQH0LD3 1 0 . *A  I NT<  < QQC  l > - 1 0 . 3 /I  0.  » 

00930  980  D0  1949  J=2, 100,1 

00990  99 ( J 3 399H0LD-S  3U9H0LD3QQ( J 3 

01000  IF(.ABS(Q9(J3  -QQC  3 .LT*  ABS(l.l*S33  S30.1*S 

01010  I F( ABS( S 3 .LT.  0.0001  3 G0  T0  1050 

01029  I F( QQC  .LT.  0.3  G0  T0  1040 

01030  I F ( QO ( J > .GT.  Q9MAX3  G0  T0  1050 
01040  1040  I F( Q9( J3  .LT.  QQMIN3  G0  T0  1050 
01050  1050  00  1330  K3!, 100,1 

01060  I F(  N .EQ.  13  QQSA  VE(  K 3 3QQ(K3 

01070  I F( QQ( K 3 .LT.  QQC  .AND.  N .EQ.  23  G0  T0  1330 
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01080  IF<*  .EQ.  u+\*i  G0  T0  13  40 
01090  Q=PMQ'J<K)/1>30> 

01100  I FCAOS < QQ <K  > > .GT.  0.001)  G0  T0  1110  I RX0VKY a 1 00000 . » G0  T0  1120 
01  l 10  1110  K X 0 VK  Y 3 ( ( 1 • ♦ Ml  A*C0S<Q)+M1  T+SIMCQ) ) /(KRA+SI  N<  3)  > -M2T)  AM2A 
01120  1120  ORsATAMtC  1 . -'ll  A *M1  A)  *SI  NO) /<  < 1 . +M  l A*M1  A)  *C0SO>  ♦ ? • *M1  A)  > 
01130  G=KRA*SI  'l(0)/tl  . ♦,'11  A + C0SO)  *M l T* S 1 MC Q ) ) 

01140  A»tl  . -G*M2T)*C1  . -G*.M2T)-G+G  I B3 -2 . *G*M2A*C  l . -G+M2  T) 

01 1 50  C3-G*G  + < 1 . -M2A*M2A>  I 003< -B-S1GM< 1 . , B)*SQRT< 8*B- 4. *A*C ) ) /2 . 

01160  QT3ATAN2(00»A)  I OQT3 1 80 . *0 T/P I 

01170  T1 =ARA*(M2A+C0S(Q D > J VI =ARA  I R1 »C05< QR> -Ml  A »C1 *M1 A+C0S< 9> 

01180  T2»RH3RA*<  1 .♦M2*C0S<  BETA-QT)  ) » V2=RH0RA*(M2A-KX0VKY*M2 T) 

01190  R23M1*C0S< ALPWA+QR) -1 . )C23! • ♦Ml +C0S( ALPHA-Q) 

01200  T33SIM(QT-RETA>  f V3  3 -<  S I N<  BETA  > + XX0  VKY*C0S(  BE  TA  ) > )R3*0.»C330. 
01210  0=*  T!  +V2*R3+V1  *R2*T3+R!  * T2*  V3-R1  *V2*T3-T1  *R2*V3-V1  *T2*R3 
01220  T«<C1 *V2*R3^V1 *R2*C3+R1  *G2* V3-R1 * V2+C3-C 1 *R2* V3- VI *C2* R3) /D 
01230  V3<  T1  +C2*R3«,C1  *R2*T3+RI  +T2*G 3-R1  *08* T3 - T1  +R2+C3-C  1 * T2*R3 > /O 
01240  R3<  Ti *V2*C  3*V1 *C2 * T3  + C 1 * T2* V3 -Cl  * V2* T3- T1 *C2* V3- VI  ♦ T2tR3 ) /D 
01250  VI A=M1 A*A1  )V1T  = M1T*A1  JV2A=M2A*A2  % V2T*M2 T*A2 
01260  Z1  »(AI  + V1  A*C0SO)  + V1  T+S I N(Q ) ) *< A 1*C 0S( 9 ) ♦ VLA) 

01270  Z2  = CA2+ V2A*C0S(QT)  + V2T*SIN(  ,QT)  )*(  A2*C0S(QT)  + V2A) 

01260  ATTEN=-4.342945*AL0GC  T* T*RH0RA*ARA*ARA*ARA*ABS( Z2/Z1 ) ) 

01290  ATTEN2<K)=ATTEN2(K)+4TTEN 
01300  G0  T0  < 1 310. 1320), M 

01310  1310  PRINT  1450.0a<K).Q0T.T.V. R»  ATTEN  i G0  T0  1330 

01320  1320  PRIMT  1440.  QQ< K)  . OQT. T. V* R» ATTEN, A TTEN2< K ) » 9QSA VE< X ) 

01330  1330  QQ(K)=QQT 

01340  1340  G0  Tq  C 1350.  1400).  M 

01350  1350  I F< INPUT  .EQ.  2)  G0  T0  1380 

01360  VI A = V2A  J V2A  = V3A  IPSI1-PSI2  JPS123PS13 

01370  PS1=PS2  J PS2aPS3  )TZ13TZ2  I TZ23TZ3  IIJW»UW23  » G0  T0  440 
01380  1380  MIA»M3A  JM2A=M4A  IM1T*M3T  )M2T»M4T 
01390  TEMP  1 * TEMP  3 J TE.MP2*TEMP4  IPS13PS3  1PS2-PS4  I G0  T0  440 
01400  1 400  PRIMT  1 460 

0141*0  PRINT:  "M0R E INPUT?  <YES*1.  N0  = 2>  " JREAD*  M0RE 

01420  G0  T0  (150,  1470),  MORE 
01430  1430  FORMA  T < FI  0 . 3 ) 

0I44P  1440  FORMA  T(  2 F9  • 3,  5F7.3,  F9.3) 

01  450  1.450  FORMA  T (2F9 . 3,  4F7.3) 

01460  1460  F0RMAT</) 

01470  1470  END 


